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O ■ Abstract 

O ■ 

For a system of n interacting electrons moving in the background of a "homoge- 
neous" potential, we show that, if the single particle Hamiltonian admits a density 
of states, so does the interacting Hamiltonian. Moreover this integrated density of 
' states coincides with that of the free electron Hamiltonian. 

c5 ■ 1 Introduction 
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We consider n interacting electrons moving in a "homogeneous" electric field in the d- 
dimensional configuration space M. d . A typical example of what we mean by a "homoge- 
neous" potential is an Anderson or alloy-type random potential. The goal of the present 
short note is to prove that, if the Hamiltonian of the single particle in the "homoge- 
neous" media admits an integrated density of states (IDS), then, so does the interacting 
n-particle Hamiltonian. Moreover, this IDS is equal to that of the interacting n-particle 
Hamiltonian. Heuristically, this is easily understood : it follows from the fact that the 
electron-electron interaction essentially lives on a strict sub-manifold of the total config- 
uration space, whereas the IDS is a thermodynamic limit over the whole space and, as 
such, measures "bulk" phenomena. 

1.1 The interacting multi-electron model 

The non-interacting n-electron Hamiltonian satisfies H n = — A +V n where the Laplacian 
— A on M. nd describes the free kinetic energy of the n electrons. As all the electrons are 
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in the same background, the potential V n is of the form 

n 

fc=l 

Hence, the noninteracting n-electron Hamiltonian takes the form 

n 

H n := l'Sfetf ® #i ® I^m"^ where i?i = -A + V x . (1.2) 

it=i 

On the one particle potential V\, we assume that 

(H.l.a) (Vi) + := max{Vi, 0} is locally square integrable and (Vi)_ := max{— V\, 0} is an 
infinitesimally — A-bounded potential that is X>((Vi)_) D T>{— A) and for all a > 0, 
there exists 7(a) < 00, such that for all G T>{— A) 

IKKMH < Qr|| A ^||+ 7 (a) W (1.3) 

(H.l.b) the operator ifi admits an integrated density of states, say Ni, that is, if i?i|A(o,L) 
denotes the Dirichlet restriction of iJi to a cube A(0, L) centered at of side-length, 
L, then the following limit exists 

N^E) := lim L- d Trace(l ] _ oo , E] (if 1 | A(0iL) )). 

L — >+oo 

Assumption (H.l.a) implies essential self-adjointness of — A +Vi on C^°(M. d ) by [RS2] 
Theorem X.29. Indeed, 

n 

V n = V+-V-, V^(x 1 ,...,x n ):=Y^(V 1 ) ± (x j ) (1.4) 

where 

• V^" is infinitesimally —A-bounded i.e. (1.3) holds for the same constants and the 
Laplacian over WL nd , 

• is non-negative locally square integrable. 

The self-adjoint extensions of — A +Vi and — A +V n are again denoted by H 1 and H n ; 
they are bounded from below. 

Classical models for which the IDS is known to exist include periodic, quasi-periodic and 
ergodic random Schrodinger operators (see e.g. [PaFi]). 

In the definition of the density of states, we could also have considered the case of Neu- 
mann or other boundary conditions. 

The interacting n-particle Hamiltonian is of the form 

H:=-A+Vi + V n , (1.5) 

where 

Vi(x u ...,z„) := Y v i x k-x{) (1.6) 

l<k<l<n 

is a localized repulsive interaction potential generated by the electrons; so we assume 
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(H.2) V : R d -> R is measurable non-negative, locally square integrable and V tends to 
at infinity. 

The standard repulsive interaction in three dimensional space is of course the Coulomb 
interaction V(x) = l/\x\. In some cases, due to screening, it must be replaced by the 
Yukawa interaction V(x) = e~' x '/|x|. 

Finally, we make one more assumption on both V\ and V : we assume that 
(H.3) the operator Vi(H n — i) -1 is bounded. 

Assumption (H.3) is satisfied in the case of the Coulomb and Yukawa potential for those 
V\ satisfying (H.l.a): H n is self-adjoint on D(H n ) C T>(— A), hence ||^(iJ„ — i)~ 1 \\ < 
\\Vi(- A • ||(- A -%){H n - where ||(- A -%){H n - < oo due to closed 

graph theorem and ||Vi(— A — i)~ 1 \\ < oo for Coulomb and Yukawa interaction potentials 
Vi, see [RS2] Theorem X.16. 

2 The integrated density of states 

We now compute the IDS for the n-electron model. Let Ax = A(0, L) be the cube in M d 
centered at with side-length L and write A£ = x ■ ■ • x Al for the product of n copies 
of Ax,. We denote the restriction of the interacting n-particle Hamiltonian H to A£ with 
Dirichlet boundary conditions by H\^n. Clearly assumptions (H.2) and (H.l.a) guarantee 
that H\\n is bounded from below with compact resolvent. Hence, for any E 6 R, one 
defines the normalized counting functions 

N L {E) := L-" d Trace(l ] _ 00 , £ ;](^|A2)). 

As usual, N, the IDS of H is defined as the limit of Nl(E) when L — > +00. Equivalently, 
one can define the density of states measure applied to a test function if as the limit 
of L _nd Trace[9?(if|A£)]. If the limit exists, it defines a non-negative measure. It is a 
classical result that the existence of that limit (for all test functions) or that of Nl(E) 
are equivalent ([PaFi]). 

2.1 The IDS for the noninteracting n-electron system 

Recall that, by assumption (H.l.b), the single particle model Hi admits an IDS (see [PaFi]) 

and a density of states measure denoted respectively by Ni and V\. 

Let H n \\™ be the restriction of H n to A£ with Dirichlet boundary conditions. One has 

Lemma 2.1. The IDS for the noninteracting n-electron model given by 

N n {E) := hm -^Trace(l ] _ 00 , £]j ff n | AZ ) (2.1) 

exists and satisfies 

N n = Nx * 2/1 * • • • * ul (2.2) 
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Let us comment on this result. First, the convolution product in (2.2) makes sense as all 
the measures and functions are supported on half-axes of the form [a, +00); this results 
from assumption (H.l.a) . When the field V\ is not bounded from below, one will need 
some estimate on the decay of N\ and v\ near —00 to make sense of (2.2) (and to prove 
it); such estimates are known for some models (see e.g. [PaFi, KlPa]). 

Proof of Lemma 2.1. The operator H n is the sum of n commuting Hamiltonians, each 
of which is unitarily equivalent to Hi, so is i? n |A n , its restriction to the cube A£. As the 
sum decomposition of H n commutes with the restriction to A£, the eigenvalues of H n \\™ 
are exactly the sum of n eigenvalues of Hi restricted to A^. This immediately yields that 

Trace(l ] _ 00 , E] (ff n | A »)) = (A>f * v[ * ■ ■ ■ * %){E) 

where N[(E) is the eigenvalue counting function for Hi restricted to Al, and v\ its 
counting measure (i.e. dN±). The normalized counting function and measure, N[ and 
v\ , are defined as 

Nf = ±fi} and ^ = ^A- 

The existence of the density of states of Hi then exactly says that N[ and v\ converge 
respectively to Ni and v\. The convergence of * v\ * ■ ■ ■ * v\ to N\ * v\ * ■ ■ ■ * v\ is then 
guaranteed as the convolution is bilinear bi-continuous operation on distributions. This 
completes the proof of Lemma 2.1. □ 

Let us now say a word on the boundary conditions chosen to define the IDS. Here, 
we chose to define it as a thermodynamic limit of the normalized counting for Dirichlet 
eigenvalues. Clearly, if we know that the single particle Hamiltonian has a IDS defined 
as the thermodynamic limit of the normalized counting for Neumann eigenvalues, so 
does the non-interacting n-body Hamiltonian. Moreover, in the case when the two limits 
coincide for the one-body Hamiltonian, they also coincide for the non-interacting n-body 
Hamiltonian. Using Dirichlet-Neumann bracketing, one then sees that the integrated 
densities of states for both the one-body and non-interacting n-body Hamiltonian for 
positive mixed boundary conditions also exist and coincide with that defined with either 
Dirichlet or Neumann boundary conditions. 

2.2 Existence of the IDS for the interacting n-electron system 

Let H\\n denote the restriction of H to the box A£ with Dirichlet boundary conditions. 
Our main result is 

Theorem 2.2. Assume (H.l), (H.2) and (H.3) are satisfied. For any E Cg°(R), one 
has 

±-J iace [cp(H lAl ) - cp(H nlAl )} 0. (2.3) 
As the density of states measure of H is defined by 



(p,dN) = ^hm^-^Trace^^iAn)]. 



we immediately get the 
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Corollary 2.3. Assume (H.l), (H.2) and (H.3) are satisfied. The IDS for the interacting 
n- electron model H exists and coincides with that of the noninteracting model H ; hence, 
it satisfies 

N = N n = Ni * v\ * ■ ■ • * v\. 

Note that, in view of the remark concluding section 2.1, we see that the integrated 
density of states of the interacting n-body Hamiltonian is independent of the boundary 
conditions if that of the one-body Hamiltonian is. 

One of the interesting properties of the integrated density of states is its regularity; 
it is well known to play an important role in the theory of localization for random one- 
particle models (see e.g. [Stl]). Usually, it comes into play through a Wegner estimate i.e. 
an estimate of the type 

E(Tracel ]£ ; _ r;ijE ; 0+r? [(if| A n)) < C w v\K n L \ 

For a specific model of random one-particle Hamiltonian, a Wegner estimate for the IDS 
of the interacting Hamiltonian was proved in [Zc]. This estimate yields Lipshitz continuity 
of the IDS for that model. 

On the other hand, Corollary 2.3 directly relates the regularity of the IDS of the 
interacting system to that of the IDS of the single particle Hamiltonian. The regularity 
of the IDS of the single particle has been the subject of a lot of interest recently (see 
e.g. [CHK, Stz]). 

Proof of Theorem 2.2. We take some q > ^ and specify the appropriate choice later on. 
By assumption (H.l. a) and (H.2), there exists ( > such that 

-oo < -C < min (inf {inf [a{H n \ Al ) U a{H\ A n)]}, mi[a(H n ) U a{H)]j . (2.4) 

Let 7 = 7(1/2) be given by (1.3) for a = 1/2. Fix A > C + 27+ 1. 

By (2.4), we only need to prove (2.3) for <p e Cq°(M) supported in (— £— 1, +oo). For such a 
function, let <p be an almost analytic extension of the function x i— > (x+\ ) q ip(x) G C£°(R) 
i.e. Cp satisfies 

• ^ e S({z e C : \%z\ < 1} 

• for any k G N, the family of functions (x i— > + iy)\y\~ k )o < \ y \ < i is bounded in 
S(R). 

The functional calculus based on the Helffer-Sjostrand formula implies 

<p{H\ A n) - Lp(H n \ A n) = (2.5) 

(H n]A n + \ )- q (H n]A n - z)- l ]dz A dz. 

In the following, we apply an idea, which has already been used in [Kl] and [KlPa] and 
which simplifies in this situation. Using resolvent equality, the integrand in (2.5) is written 
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as 



(H {A n + Ao)- 9 (if|A2 - zY 1 - (H n[Al + A )^ 9 (if„|A2 - z)' 1 = 

= (H n{Al + \ y q [(H lAl - z y l - (H n \ Al - z)- 1 ] 

+ [{H\ A n + A r f/ - (H n]A n + X )- q ](H ]A n - Z )~ l 
= -{H n \ A n + \ Q )-1(H n \ Al - Z^ViiH^n - ZY 1 



— / j(H n \ A n + Ac))' q l Vi(H\ A n + A ) \H\ A n 
1=1 

(2.6) 

Estimating the trace of (2.6), we choose e > and write 

V i = V l - l m < e] + Vi ■ l mi >e } (2.7) 

and note, that Vi ■ l{|y i i< £ } is bounded by \\Vi ■ l{|y i |< £ }|| < £■ As V is non-negative, one 
has 

Slippy ■ l { |v < |> e} ) C |J |J J (x 1; . . . ,x n ) G M nd : V(x t - Xj ) > £ _ I . (2.8) 

As, by assumption (H.2), V tends to at infinity, (2.8) implies that there exists < C(n; e) 
(independent of L) such that 

n{{\Vi\ >£}nA2) <C(n,e)L^ d , (2.9) 
where /x(-) denotes the Lebesgue measure. Using decomposition (2.7) of Vi, we obtain 
|Trace(i7 n | A « + X Q )- q (H nlA n - zyW^H^n - z)" 1 ! 
< -^-Trace|(i/ n | A n+Ao)^| 



'\Vi(H\ A n - z) Trace] (H n \ A n + A ) 9 l{|V;|> £ }nA£| 



\$sz 



• ||(#n|A£ + Ao)~ 1 l{H/,|> £ }nA2llr £( • HV^ifnlA" + A ) _1 || 

(2.10) 

where || • ||r 9 denotes the g-th Schatten class norm (see [Si]) and we used Holder's inequality. 
In the same way, the cyclicity of the trace yields 

|Trace(# n | A « L + X^-^V^H^ + X )~\H\ Al - z)' 1 ] 

< Trace|(if| A2 + X o y l (H n]A n + \ ) l -«- l V t {H\ Al - z)' 1 ] 

< \\(H ]A n + X Q )- l (H n]A n + A ) l || ■ Trace\(H n \ A n + Ao)^ 1 ^^ - z)~ l \ 



Q i 
< 7^— r||(-Hn|A2 + \))~ ||r ' ll(-^n|A£ + -V>)~ l{|Vi|>e}nA£ \W q ' ||Vi(-ff n |A£ + ^o)~ I 



+ ^ll(^|A2 + A o)- 1 ||^ 



(2.11) 



We are now left with estimating || (H n \ A n + A ) 1 \\ Tq and || (H n \ A n + A ) 1 l { \ Vi \ >£}nA i \\ Tq for 
q sufficiently large, depending on nd. Therefore, we compute 



(-^n|A™ + Ao) 1 l{|y i |> e }nA™||r 9 < \\{H n \ A -n + A ) 1 (— A| A n +A ) 2 \\r 2q 

■ || (- A|A 2 +Ao)~2l { | yi | >e}nA n|| T2(j (2.12) 



We use the decomposition (1.4). As the Laplacians are positive, the infinitesimal — A- 
boundedness on V~, Theorem X.18 of [RS2] and the definition of 7 imply the following 
form bound 

K0,^a^>I<^-A|a 2 <P)+lU\\ 2 - 

As A > 27 + 1, one has 

H n \ Al + A > - A| A » +V~ lAl + A > -(- A| A » -2 7 + 2A ) > -(- A, A n +A ). 
Thus, the operator H n \ A n 4- Ao is invertible and 

{H n \ A n + A ) _1 < 2(— A| A n +A ) _1 . 

Let and (<j>j)j respectively denote the eigenvalues and eigenf unctions of the Dirichlet 
Laplacian — A| A « (the index j runs over (N nd )*). For q G N such that 2q > nd, we 
compute 

||(# n | A « + Ao)-^- A| A n +A )5||| 9 = M" A |A£) + Ao) 9 %, (H n]A1 + Ao)" 1 ^) 29 

< 2 2 * J2 M- A \Al) + X o) q (^ (" A, A n +Ao)^V,)) 29 

= 2 2 « Yl M~ A \Al) + Ao)" 9 < CL nd . 

The last estimate is a direct computation using the explicit form of the Dirichlet eigen- 
values. 

By Lemma 2.2 in [KlPa], we know that, for q e N such that 2q > nd, there exists C q > 
such that, for any measurable subset A' C A£, one has 

||(- A| A n +X )-h A .\\% t < C q a{M). (2.13) 

Choosing A' = {| Vi\ > e} H A£ and taking (2.9) into account, then by combining estimates 
(2.10)-(2.13), we get that there exists c, depending only on q (and the bound in assumption 
(H.3)), such that 

Trace|(if| AZ + X y q (H ]Al - z)~ l - (H n]A n + \ )-«(H n]A1 - z)- l \ 

\<sz\ z \ksz\ z \<sz\ \ksz\ 

Using this inequality in (2.5), we get (2.3) as, (p being almost analytic, d(p{z) vanishes 
to any order in Qz as z approaches the real line. Thus, we completed the proof of 
Theorem 2.2. □ 
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